Singular surfaces and cusps in symmetric planar 3-RPR manipulators 
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Abstract — We study in this paper a class of 3-RPR ma- 
nipulators for which the direct Jiinematic problem (DKP) is 
split into a cubic problem followed by a quadratic one. These 
manipulators are geometrically characterized by the fact that 
the moving triangle is the image of the base triangle by an 
indirect isometry. We introduce a specific coordinate system 
adapted to this geometric feature and which is also well adapted 
to the splitting of the DKP. This allows us to obtain easily precise 
descriptions of the singularities and of the cusp edges. These 
latter second order singularities are important for nonsingular 
assembly mode changing. We show how to sort assembly modes 
and use this sorting for motion planning in the joint space. 

I. Introduction 

Planar parallel manipulators have received a lot of atten- 
tion [1], [2], [4]-[19] because of their relative simplicity 
with respect to their spatial counterparts. Moreover, studying 
the former may help understand the latter. Planar manipu- 
lators with three extensible leg rods, referred to as 3-RPR 
manipulators, have often been studied. Such manipulators 
may have up to six assembly modes (AM) [2] and their 
direct kinematics can be written in a polynomial of degree 
six [3]. It was first pointed out that to move from one 
assembly mode to another, the manipulator should cross 
a singularity [2]. However, [4] showed, using numerical 
experiments, that this statement is not true in general. More 
precisely, this statement is only true under some special 
geometric conditions, such as similar base and mobile plat- 
forms [5], [6]. Recently, [7] provided a mathematical proof 
of the decomposition of the workspace into two aspects 
(singularity-free regions) using geometric properties of the 
singularity surfaces. Since a parallel manipulator becomes 
uncontrollable on a singular configuration, the possibility to 
change its assembly-mode without encountering a singularity 
is interesting as it can enlarge its usable workspace. Knowing 
whether a parallel manipulator has this feature is of interest 
for both the designer and the end-user. The second-order 
singularities, which form cusp points in plane sections of the 
joint space, play an important role in non-singular assembly- 
mode changing motions. Indeed, encircling a cusp point 
makes it possible to execute such motions [5], [9], [12]- 
[19] A special class of planar 3-RPR manipulators has 
been studied recently [9], [10]. These manipulators have 
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the peculiarity that the resolution of the direct kinematics 
problem is split into a cubic equation and a quadratic 
equation. Their geometry is characterized by the fact that 
the platform triangle is congruent to the base triangle via an 
indirect isometry of the plane; this is the reason why we call 
them "symmetric". 

We propose here a coordinate system for the workspace 
which is adapted to this specific class and reflects the 
splitting of the direct kinematic problem (section II). We pay 
attention to the description of singularities (section III) and 
cusps (section IV) using these coordinates. We show how 
to sort assembly modes and use this sorting to do motion 
planning in the joint space (section V). 

II. Alternative coordinates for the workspace 

The base triangle is denoted by A1A2A3. In the direct 
orthonormal frame T with origin Ai and first axis oriented 
by A1A2, the coordinates of A2 are (&, 0) and those of A3 
are [d, h). The platform triangle is denoted by B1B2B3. Due 
to the symmetry property, the coordinates of B2 and B^ in 
the direct or thono rmal frame with origin Bi and first axis 
oriented by B1B2 are respectively (&, 0) and {d,—h). The 
length of the leg AjBi is as usual denoted by pi. 




Fig. 1. A symmetric 3-RPR manipulator 

The platform triangle B1B2B3 is the image of the base 
triangle A1A2A3 by a glide reflection S. We encode this 
glide reflection S by the triple {il),r,g) such that the glide 
reflection is the orthogonal symmetry with respect to the 
line A with equation x cos(V') + y sin(-0) — r = followed 

— sin('0) 

cos(-!/') 

symmetry axis (the equation of A and the coordinates of the 
translation vector are given in frame F attached to the base 
triangle - see figure |2]i. 

We choose the angle ip in [—it/2, 7r/2] and make the 
identification of (— 7r/2, r, g) with (7r/2, — r, — g). 



by the translation of vector 2g 



parallel to the 




Fig. 2. The glide reflection sending AiA2A's to B1B2B3. 



Usually the workspace is viewed as the space of rigid 
motions in the plane and a pose of the manipulator is encoded 
by the rigid motion R carrying the half-line [^1^2) to the 
half-line [BiB2)- The rigid motion R and the glide reflection 
S are related in the following way: 5* is the orthogonal 
symmetry with respect to {A1A2) followed by R. If the rigid 
motion R is given (as in [9], for instance) by the angle of 

rotation ip and the translation vector ( 1 , then the relation 

\ y J 

between the two systems of coordinates is as follows: 

(f = 2i/j + Tr (mod 2tt) (1) 

X ~ 2 (r cos('i/;) — g sm{tp)) (2) 

y = 2(r sin(i/')+5 cos(i/')) • (3) 

It is easy to compute the lengths pi of the legs AiBi in 
terms of {ijj,r,g), since Bi is the image of Ai by the glide 
reflection. The square pf is the sum of the square of the 
double of the distance of Ai to the axis A and the square of 
the norm of the translation vector, which is 4g^. This gives: 



pI 
pI 
pI 



- 4(r2+.g2) , (4) 

= 4 ((6 cos(V;) - rf + g") , (5) 

= 4 ((d cos(?A) + h sin(^) - rf + g^) . (6) 

It will be convenient to introduce 82 =■ {p\ ~ Pi)/4 and 
(^3 = (pg — p?)/4. These quantities depend only on •(/; and r, 
and not on g: 

82 = -2 6 cos (V') r + 52 (cos {il^)f , (7) 

5s, = (dcos(V') + /isin(V'))^-2r(dcos(-0) + /isin(V')X8) 

Eliminating r between these two equations and writing the 
equation obtained in i = tan(V'), we get the third degree 
equation: 



52ht^ + {bh^ - b53 + S2d) <2+ 
(2 bdh - b^h + S2h) t-b53 + S2d 



bd^ 



b^d 



0. (9) 



We shall use {p'i,P2,P3) as coordinates for the actuated 
joint space. Of course, this joint space is contained in the 
positive orthant (IR+)'^. 

The direct kinematic problem (DKP) can be solved as 
follows: 

1) Take a real solution in t of equation (|9]l, which deter- 
mines ip = arctan(i). Generically there are 3 or 1 real 
solutions, depending on the sign of the discriminant of 
the equation. 

2) Compute r from ip using equation (|7]), which gives 



b cos (■(/)) 



(10) 



b cos(-0) ^ 

3) Solve the equation (|4|i for g. It has two real opposite 
solutions when p^ > Ar^. 

One should take care in the resolution of the DKP of 
the case -0 = ±7r/2, i.e. t = 00. This corresponds to the 
vanishing of the third degree term in equation ^, which 
occurs when 82 = 0. In this case one can use equation ^ 
to compute r, which gives r = {h? — S3)/2h for ij.' = tt/2. 

Note that the existence of a solution to the DKP depends 

> first, on the existence of a solution (-0, r) to the system 
of equations (7) and (8), 

> second, given such a solution {ip,r), on the existence 
of a solution g to equation (4). 

These two conditions of existence will be important for the 
discussion of singularities in the next section. 

III. Singularities 

The singular surface in the actuated joint space is thus 
given as the union of two surfaces 5*1 and ^2, corresponding 
respectively to steps 1 and 3 of the resolution of the DKP 
described above. The fact that the singular surface splits in 
two components has already been observed in [9]. We will 
now describe these two surfaces. We will also describe the 
critical surfaces Si and S2 in the workspace, whose images 
by the mapping {ip,r,g) n> {pi,P2,P3) given by equations 
dUi, Q and (|6]l are ^i and ^2 respectively. 

A. The first singular surface 

The surface Si is the intersection of the actuated joint 
space (always with coordinates {p^, P2i P3)) with a cylinder 
having generatrix parallel to (1,1,1) and basis a curve C in 
the plane of coordinates {52,53). An equation for C can be 
obtained as the discriminant of equation (|9|i; it is a quartic. 
An alternative way to describe C is to compute the jacobian 
determinant of the mapping $ : (0, r) M> (52,(53) given by 
equations ^ and ^. The jacobian curve F in the space with 
coordinates {ip,r) is given by 



cos(0) 



(h^ + bd- d^) cos(V') sin(V') + 
(2d -b)h cositPf + {b-d)h 



(11) 



This equation is essentially the same as the third degree 
characteristic polynomial obtained in [10]. 



Observe that r(?/' + tt) = —r{ip). The critical surface Ei 
in the workspace is the set of all {ip,r,g) such that (^jt) 
belongs to F. The curve C is the image of F by the mapping 



$, and it can be parameterized by rational functions of t 

tan(-0) as 



S. = 



b {{2d - b)ht^ + {2<f - 2bd - 2h^) t+{b- 2d)h) 



{ht + df{he 



{l + f^Yh 
2{d-b)t- h) 



(1 + ^2)2/1 



(12) 



So C is indeed a rational quartic. Its singular points are three 
real cusps that can be found by looking at the stationary 
points of the parameterization. These stationary points cor- 
respond to parameters t which are roots of the cubic equation 

{b-2d)ht^+'3{h^-d^+db)t^+3h{2d-b)t+d^-db-h^ = 

(13) 
Since the discriminant 108 (d'^ + fi^)'^ {{d-b^ + h^Y of this 
cubic equation is strictly positive, there are always three real 
roots and hence three real cusps on the curve C. Actually, 
transforming the equation to an equation in ip, one obtains 



tan(37/') 



d^ -bd- h^ 
{b ~ 2d) h 



(14) 



Note that the curve C has no other singular point than the 
three cusps. Indeed, a rational quartic may have only up to 
three singular points. The curve C always has the shape of 
a deltoid, i.e. a closed curve with three cusps connected by 
arcs concave to the exterior. 

B. The second singular surface 

The second critical surface S2 in the workspace is simply 
given hy g = 0; so it is independent of the geometry of the 
manipulator (this is already observed in [9]). Its image ^2 
in the actuated joint space is parameterized by substituting 
g = in equations (4-6). So the surface ^2 is the image of 
the elliptic cylinder 



Pi = 


= 2r, 




(15) 


P2 = 


= 2 (r — 6 cos(V')) 




(16) 


P3 = 


= 2{r-d cos(V') - 


- /isin(-0)) 


(17) 



by the mapping (pi,p2,/53) ^ {Pi,pI-,pI)- The implicit 
equation of ^2 can also be obtained by eliminating t between 
equation (|9]l and the equation 



Pi — 



62 



2^2 



l+t2 



where the right hand side is the expression for 4r2 derived 
from ( [Tol l. The implicit equation for 5*2 obtained in this way 
is a quartic equation in pi, P2, P3, not a very nice one. 

C. An example 

We consider the manipulator with parameters & = 1, h = 
1, d = 0. In this case the curve C in the plane {621 5^) is 
a hypocycloid with three cusps (a deltoid) inscribed in the 
circle with center (1/4, 1/4) and radius -^/O/S. The three 

cusps on C correspond to the values —--7-, 7:, -7 of ip. 
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Fig. 3. The curve C 

The critical surface S2 in the workspace is always given 
by g = 0. The critical surface Si is parameterized by 

ip, — cos(7/') (— 1 — sin(^) cos(?/') + cos(?/')2), gj . 



g 
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Fig. 4. The critical surfaces Si (in blue) and E2 (in red). 

Both critical surfaces are represented in figure H] It may be 
interesting to compare this figure with figure 2 in [9], which 
represents the same surfaces (with the same color code), but 
in a different coordinate system. The choice of coordinates 
made here "straightens" the critical surfaces. 

The three black lines of figure |4] are the lines of points 
which correspond to cusps in the joint space. 

We represent now the singular surfaces 5*1 and S2 in the 
joint space (See figures |5] and |6] ). The surface Si is a part 
of a cylinder on the hypocycloid and has three half-lines of 
cusps. The drawing of the singular surfaces is made using 
their parameterizations by {ip,g) for Si and by {ip,r) for 
S2. 





Fig. 5. Si (in blue) and 52 (in red) cut at pj = 1. 





the number of cusps is the number of cusp points in the slice 
of the joint space. Over each of these cusps there are two 
triple solutions of the DKP, corresponding to opposite values 
of g. 



P, 




Fig. 6. Si (in blue) and 52 (in red) cut at pj = 4. 



(a) 



IV. Cusps 

The singular surface Si has three half-lines of cusps, all 
parallel to the vector (1,1,1). So the cusps are entirely 
determined once we know the origins of these three half- 
lines. The three possible angles ip ^e determined by equation 
dUli: 



V'cusp = -r arctan 



(f -bd- h"^ 
{b - 2d) h 



+ k — (mod tt) 



for k — 0,1,2. For each of these three values of -iAcusp, 
the corresponding value Tcusp is given by equation (fTTI) . The 
couple (V'cusp, ''cusp) determines the line supporting the half- 
line of cusps, and the origin of this half-line corresponds to 
g — 0. So we get three values for the pf of the origins 
of the three half-lines of cusps, which are the three values 
for 4:r^ . These three values for pf are the roots of a third 
degree polynomial with coefficients depending on 6, h, d; the 
constant term of this polynomial is 

4h^ {^2d+bf {d^-2bd^+h^d+b^d-2bh^f {-h^+2bd-d^f 

and its discriminant is always nonnegative. 

Let < /?! < /32 < /33 be the three bifurcation values of 
p\ for the number of cusps, i.e. the p\ of the origins of the 
half-lines of cusps. Then the slice at p^ = c has cusp if 
< c < /?i, 1 if /3i < c < /32, 2 if /32 < c < /?3 and 3 
if Pz < c. One of the bounded intervals may be empty, if 
the constant term of the equation of the third degree in pf 
vanishes (for the first interval) or if its discriminant vanishes 
(for the second or third interval). One has to understand that 




(b) 



' Fig. 7. The bifurcation values of p^ for the number of cusps. 

Figure |7^ represents the bifurcation values as function of 
h and d, with b fixed equal to 1 (green for /3i, yellow for 
/?2 and red for /S^). Figure |7}5 shows detail for the first and 
second bifurcation values when < d < 1. 

V. Sorting assembly modes and motion planning 
IN the joint space 

The essential idea here is the following: when one starts 
from a nonsingular solution of the DKP at a point in the 



joint space with coordinates (pij/O^jPi) and moves in the 
direction of the vector (1,1,1), then the solution of the 
DKP follows smoothly, without crossing a singularity in 
the workspace. Indeed, consider equations (4-6): the motion 
pi '~^ Pi ~^ ■^^^ increasing A, can be lifted in the workspace 
(with coordinates {ip,r,g)) to a path with ijj and r fixed and 
g increasing to g^ + A^/4. 

The segment in the joint space can cross the second 
singular surface 5*2. This corresponds to the appearance of 
two new solutions to the DKP (one in each aspect), with a 
different couple {ip,r). But it never crosses the first singular 
surface S'l which is a cylinder with generatrix parallel to 
(1,1,1). 

We denote pi + P2 + Pa by i^. Then (1^,(52,(53) form a 
system of coordinates for the joint space which is convenient 
for our present discussion. Moving in the joint space in the 
direction of the vector (1,1,1) is increasing v, keeping 62 
and 63 fixed. 




-1 -0.5 0.5 1 1.5 

Fig. 8. The section u = 8. 



The situation in the section v = k stabilizes for sufficiently 
large values of k: the section of S2 is a big oval surrounding 
the section of 5*1 which is a deltoid with three cusps (the 
curve C, base of the cylinder). Inside this curve C there are, 
in each aspect, three continuous solutions of the DKP and 
between this curve C and the section of S2 there is in the 
same aspect one continuous solution of the DKP. Label by 
this solution, and label by 1,2,3 the three arcs of C between 
the cusps. Then we can label the three solutions inside the 
deltoid C by 1, 2, 3 according to the label of the arc of the 
deltoid through which they are connected with the solution 
0. Figure [8] illustrates the labelling in the same example as 
above (the coordinates ((52, i53) are used in the section v = 8). 

In this way we can label every solution of the DKP 
contained in one aspect by one of the labels 0,1,2 or 3. In 
each aspect, all points in the same label form a connected 
region and the boundaries between these regions are the so- 
called "characteristic surfaces" obtained by pulling back the 
singular surface Si in the aspect [20]. The characteristic 
surfaces in the workspace with coordinates {^',r,g) are 
cylinders with generatrix parallel to the g-axis and basis the 



two curves in the {ip, r) plane which are given by 

h{2 d + b) cos {tp) + {K^ + bd- d^) sin (V) 



Ah 



± 



^/{h^llFY(h?Td?T¥~2bd) 

Ih 



Figure |9] represents a section of the characteristic surfaces 
(in green) by a plane g ~ constant of the workspace. The 
blue curve is a section of Si and separates the aspects. In 
each aspect, the characteristic surface delimitates the four 
sorts of assembly modes. The label correspond to points 
which are mapped outside of the deltoid C in the joint space, 
and the labels 1,2,3 to points which are mapped inside. A 
path from to 1 inside an aspect is mapped to a path going 
through the arc of the deltoid with label 1, etc.. 




Fig. 9. The four sorts of assembly modes in each aspect, in a section g = 
constant of the workspace. The green curves are sections of the characteristic 
surfaces. 



We illustrate how the labelling can be used for motion 
planning in the joint space with an example, again for the 
manipulator with parameters b = I, h = 1, d = 0. We choose 
a goal position for the manipulator, given by ?/> = 7r/4, r = 
1.1 and 3 = 0.4 (Figure [JOll. 
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Fig. 10. The goal position. 

The goal position corresponds to values pi k, 2.34, 
/52 = Pa ~ 1-12 and v k, 7.99; it is labelled and mapped 
outside of the deltoid in the joint space. We explain how 
to plan (in the joint space) a path to the goal position 



from a position of the manipulator in the same aspect, with 
V < 7.99. The starting position corresponds to p", pjj P3 and 
label £e {0,1,2,3}. 

> Increase simultaneously pi,p2, ps following pi = 

V(pl')' + t2until;. = EP?=8 
« Keeping ly = 8 constant, move in the plane with 
coordinates {62,53) from 

{s'2 = ip'^r (p;)^ ^i = iptr (p?)^) 

to (—1.06, —1.06) following a path inside the red curve 
and crossing only arc ^£ of the deltoid if the label is 

e. 

Figure [TT] shows such a path for label 1 = 3. (Actually, 
it shows only the part of the path in the plane (62,63) for 
1/ = 8, since the first segment of the path increases ly without 
changing 62 nor 63). 




-1 -0.5 0.5 1 1.5 

Fig. 11. A path from a starting and a goal position 



VI. Conclusions 

The choice of coordinates for the workspace well adapted 
to the special class of symmetric manipulators allowed us 
to take full advantage of the de-coupling of DKP (into a 
cubic and a quadratic equation) in the computations. We 
obtained rational parameterizations of the singular surfaces 
in the joint space. We obtained also a good description of the 
cusp curves on these surfaces, as well as precise information 
on the bifurcation of the number of cusps in the slices of the 
joint space. Finally, we were able to sort assembly modes in 
an aspect and to use this sorting for motion planning in the 
joint space. 
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SOO based algorithm 1 



Get current posture q 

Initialize the trajectory Traj by adding q to tail 

Current position of the end effector P_end=direct_kinematic(q ) 

Distance D=Dist(P_end, P_des) 

while D>/j 

1 . Generate a set of new postures in 

the neighbourhood (q±dq) 

2 .Select the collision-free postures 
Calculate the distances Dist_set 

3. Find out the nearest collision-free posture 
q=min(Dist_set) 

4. Update current posture q 

5. Update the distance D 

6. Add the cuirent new posture q to trajectory Traj 
end while 



i* 
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